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I revue the so called Wilson loop approach to bound state problem in QCD. I shall 
show how using appropriate path integral representations for the quark propaga- 
tor in an external field it is possible to obtain corresponding path integral rep- 
resentations for various types of gauge invariant Green functions which have the 
important feature of involving the gauge field only trough Wilson loop correlators 
or their generalizations. Two different kinds of representations are used, one given 
in the form of a semi-relativistic expansion, the second completely relativistic of 
the Feynmann-Schwinger type. In this way starting from reasonable ansatz on the 
non perturbative part of the Wilson correlator one can obtain: expressions for the 
semi relativistic (spin dependent and momentum dependent) qq and 3q potentials, 
a "second order" qq Bethe-Salpeter equation and and a related Dyson-Schwinger 
equation. I shall concentrate on the three quark potential for which new contro- 
versial results have been obtained by lattice numerical simulations and on a three 
dimensional reduction of the BS equation obtained in the form of the eigenvalue 
equation of of a squared or a usual mass operator. We shall report on a numerical 
resolution of such equations which seems to give a comprehensive reproduction of 
the entire meson spectrum with the exception of light pseudo-scalar bound states 
for which a complete four dimensional treatment should be necessary. 



1 Green Functions 

The single quark, the quark- antiquark and the three quark gauge invariant 
Green functions can be written as 

G gi (x -y) = (0|T£%, x)V(a#(2/)^(y)|0) = Tv c (U(y, x)S(x, y; A)) , (1) 

G^(x 1 ,x 2 ;y 1 ,y 2 ) = ^(0\T^ 2 (x 2 )U(x 2 , xM^i(yiP(yu ifc)^(ife)|0) - 

= ^Tr c (U{x 2 ,x 1 )S 1 (x 1 ,y 1 ;A)U(y 1 ,y 2 )S 2 (y 2 ,x 2 ;~A)) (2) 

(the two quarks are supposed to have a different flavor, otherwise an annihila- 
tion term should be added to the last term) 



G gl (x 1 ,x 2 ,x 3 ;yi,y 2 , y 3 ) = ■^ea 1 a a a s £b 1 b 3 ba (°I T U a3 ° 3 (x M ,x 3 )U a2C2 {x M ,x 2 ) 
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U aici (x M , zi)-03c 3 Oa)^ {x 2 )ipi ci {xi)vb ldi (yi)Tp 2d2 (y 2 )i> 3d3 (y 3 ) 
U dlb Hyi,yM)U d2b2 (y 2 ,y M )U d ^(y 3 ,yM)\0) = ^e a ia,a s £b x b 2 b 3 
({U(xM-,xi)Si (xi, J/1 ; A)U(yi,yM)) aibl (U(x M ,x 2 )S 2 (x 2 ,y 2 ;A) 
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U(y 2 , yAi)) a2b2 (U(x M ,x 3 )S 3 (x 3 ,y 3 ;A)U(y 3 ,y M )) a:>b;> ) ■ 



(3) 



To avoid complicate indices in the above equations we have identified the var- 
ious type of functions simply by their arguments; furthermore tp c denotes the 
charge-conjugate fields, and U the path-ordered gauge string (Schwinger string) 



(the integration in (Q) is along an arbitrary line joining a to 6), the tilde and Trc 
denote the transposition and the trace over the color indices alone. The second 
equalities are the result of an explicit integration over the fermionic fields and 
the angle brackets denote average on the gauge variable alone (weighted in 
principle with the determinant Mf(A) resulting from the integration). The 
quantities S, ■ ■ ■ S 3 are the quark propagators in the external gauge field A 11 , 
which are defined by equations of the type (we shall suppress indices specifying 
the quarks, as a rule, when dealing with single quark quantities) 



Here to is supposed to hold an infinitesimal negative imaginary part (to = 
to — iO). 

2 Quark-antiquark potential 

By performing a Foldy-Wouthuysen transformation on Eq. (|s|) we can replace 
the the 4x4 Dirac type matrices S(x, y; A) with a 2 x 2 Pauli propagator 
K (x, y; if , t{) satisfying a Schrodinger-like equation. Solving such equation by 
standard path-integral technique, we obtain 



Here T is the quark world line connecting (y, t\) with (x, if) over which the path 
integral acts, T s and P are time ordering prescriptions for the spin and the color 
matrices respectively. Furthermore, as usual F^ v = d^A u — d v A^ + ig[A^, A v \ , 
P^u _ i s »u P ap^ and jyp vji = Q v p vil + ig[A v ,F vli ] , e^ pa being the four- 
dimensional Ricci symbol. 




[yfDp- m)S(x,y;A) = S 4 (x - y) . 



(5) 
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Replacing @ in (|) we obtain a path integral representation for a Pauli 
type qq propagator which can be compared with the corresponding expression 
in potential theory (see P] for details). Then one find that, in order the inter- 
action between the two particles can be described in terms of a potential, a 
function V q ^(zi, z 2 , Pi, p 2 , Si, S 2 ) must exist such that op to the order 



■lnW q - q + ij2— I dx^(s l j {{F l ^x)))--^-S l j e lk ^{{F fir (x))) 
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j-j' J 

{[{Mx)F ka [x'))) - ({F ¥ (x))){{F ka (x>)))^+... = j\tV™, (7) 

with the notation ({f[A})) = (W qq f[A\) / (W qq ) and W qq is the Wilson loop 
correlator defined by 

W qq = i ^Tr c Pexp (ig jf ^(a;)^ \ . (8) 

In Eq. (||) the integration loop T qq is supposed to be made by the quark world 
line Ti, the antiquark world line T 2 described in reverse direction and the 
two Schwinger strings that close the curve and can be here taken equal time 
straight lines; the symbol P denotes color ordering along T qq . 

Notice that the quantities ((F^(zj))) and {{F lxu {zj)F prJ {z' : j))) can be ex- 
pressed as functional derivatives of i In W qq and so the potential is determined 
in principle by such quantity alone. 

Eq. (0) can be reelaborated in various way. By expanding i In W qq and 
the spin dependent terms in Zj it is possible to recast the resulting local co- 
efficients in terms of static Wilson loop correlators (straight quark world-lines 
parallel to the time axis) with insertion of field strengths. Such expressions 
are particularly suitable for numerical simulation and to this aim have been 
successfully used ft]. 

To obtain analytical expressions today we have to relay on models for 
the evaluation of the Wilson correlator. The simplest model is the so called 
minimal area law model (MAL model) and consists in writing i In W qq as the 
sum of its perturbative, an area and possibly a perimeter term 

i\nW qq = {i In W qq ) pcl . t + <rS min + \cP (9) 

where S^in denotes the minimal surface enclosed by the loop T qq and P length 
of the loop. In Eq. (H) the first quantity is supposed to give correctly the short 
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i(ln W, 5 ) port = ~-g 2 j> dz» j> dz vl D^{z - z') (10) 



range limit of the interaction and it is suggested by asymptotic freedom; the 
other two are supposed to represent the long range behavior and are suggested 
by pure lattice gauge theory and numerical simulations 

At the lowest order the perturbative term can be written as 

2 
3' 

D fll ,(z — z') being the free gauge propagator. For what concerns the area term 
it can be checked, by solving the appropriate Euler equations, that <!5 m in can be 
replaced by the surface spanned by the straight lines joining equal time points 
on Ti and T2 up to 0(v 2 ). That is, we can write 

S m in « 1' dtar f dX [1 - (Az 1T + (1 - A)z 2T ) 2 ]^ (11) 
Jti Jo 

where z,t denotes the transversal part of z_j, ij T = (5 — f h f k )z^. 

Replacing (|Io| ) and (11) in (^) and (0), by expanding in the velocities 



and other manipulations we obtain a semi-relativistic potential of the form 
V™ = K?Itic + V£ + Vfl where V* tic = -f a.1 + or + C, V s f and V* are 
certain complicate spin dependent and velocity dependent expressions which 
are reported in p] but we do not reproduce here. Similar results can be obtained 
starting from more elaborate models 



3 Three quark potential 

A semi-relativistic three quark potential can be obtained proceeding in a sim- 
ilar way on the three quark Green function (|s|). Instead of (0) this time we 
have 

i \nW 3q + ij2— I dx^(s\{{F lll {x)))-^-S l j E lkr p k j {{F lir {x)))- 



3,3' 



"3 ""3' 



((F^F^x'))) ((^(x))) ((F ka (x'))) 



dtV 3q 



(12) 



where now ((f[A])) = (Wj, q f[A\) / (W 3q ) and W 3q is the Wilson loop correlator 
for three quarks 

ai&i 



W"3</ — 2j ^aia2a3^bib 2 b 3 
a 2 b 2 



Pex.pyigJ_ dx^ 2 A fl2 (x) 



Pexp^igJ_ dx^A^ix) 
P exp (ig J_ dx^ 3 A^ 3 (x) 



asba 



(13) 



In ( |l3| ) a,j, bj are color indices, j — 1, 2, 3 and Tj denote the curve made by 
the world lines Tj for the quark j joining ijj to Xj, a straight line on the surface 
t = tj merging from an arbitrary fixed point yu and arriving to yj, another 
straight line on the surface t = tf connecting the world line to a second fixed 
x M . 

A straightforward generalization of the arguments used in the quark-antiquark 
case suggests to write in place of and ( pX| ) 

i In W 3q = | 9 2 Y, J r dx * J r dx ~J W ^( x * ~ x i) + ° S ™* + \ CP > ( 14 ) 



where the perturbative term is taken at the lowest order in a s and now S m i n 
denotes the minimum among all the surfaces made by three sheets having the 
curves Ti, T2 and T3 as contours and joining on a line Tm connecting j/m with 
xm (the minimum is understood at fixed Tj as the surfaces and Tm change ). 
Obviously, P denotes the total length of T\, T2 and T3. 

As in the qq case, up to the second order in the velocities S m i n can be 
replaced by the surface spanned by the straight lines joining the point Zm(£) 
at the time t with the equal time positions of the three quarks zi(£), ^{t) 
and Z3(i), zjvfft) being constructed according to the following rule: if no angle 
in the triangle made by Zi(t), ^{t) and ^(t) exceeds 120° (configuration I), 
Zjitf-(i) coincides with the point inside the triangle which sees the three sides 
under the same angle 120°; if one of the three angles in the triangle is > 120° 
(configuration II), Zm(4) coincides with the corresponding vertex, let us say 

The final result is of the form V 3q = + V; 3 d 9 + V; 3 d 9 , where V s 3 t q at = 

J2j<i (-§ f^) +(r(r 1 +r 2 +r 3 ) + C (with Vj = z,j-z M , r jt = r^-r; = zj -z ; ), 

while Vj, d 9 and V 3 ^ are spin dependent, and momentum depend expressions of 
order 1/m 2 for which again I refer to 

We observe also that the short range part in V 3q , coming from the first 
term in is of a pure two body type: in fact it is identical to the elec- 
tromagnetic potential among three equally charged particles but for the color 
group factor 2/3. 

The three sheet definition of S^in given above is suggested by a direct gen- 
eralization of Wilson's original argument for W q q, based on the lattice formu- 
lation. It is adopted by the majority of the authors and it is said to correspond 
to an Y flux tube configuration. However an alternative possibility has also 
been considered p] . This consists in setting 

Smin = ^Cr(5i2 + S23 + S31) , (15) 
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where now Sij denotes the minimal surface delimited I\ and Tj and the factor 
■i is introduced in order to reproduce the qq case when two quark world-lines 



being equal to the quark antiquark potential relative to the couple ij as 



collapse. Assumption J15|) is said to correspond to a A flux tube configuration. 
In this case even the long range part of the three quark potential is the sum of 
purely two body potential and one can write simply V 3q = | (V™ + V23 + V31 ) > 

defined by @. 

If we restrict ourself to the linear rising part (the confining part) of the 
static potentials alone we can write for the two configurations 

T^ onf = <T( n +r 2 +r 3 ), (16) 

VK° nt = \a{r 12 +r 23 + r 31 ). (17) 

Then, simple geometrical considerations show that 

V£ on{ < V£ onl < -^=V£ onf , (18) 
V3 

where the equal sign in the first step holds when the three quarks are aligned, in 
the second one when they form an equilateral triangle; the second one being the 
situation in which the difference between the two expressions is maximal (note, 
however that ~ 1.155). Until recently the various attempts to discriminate 
the two expressions by evaluating numerically lnW3 g for a static loop -have 
been unsuccessful. Presently preliminary results obtained by Bali et al. p] for 
the equilateral arrangement are claimed to support A configuration. On the 

fntrary the very accurate estimates reported by Matsufuru at this conference 
are definitely in favour of the Y configuration. Indeed it seems that, due 
to the large errors occurring in the results reported in Ref. p] .at the large 
distances, these can be reasonably reconciled with those of Ref. p] simply by 
an appropriate choice of the constant C (private communication by Suganuma) . 
If this were confirmed, the Y configuration should be considered established. 

4 Bethe-Salpeter and Dyson-Schwinger equation 

Full rclativistic bound state equations can be obtained along similar lines using 
a covariant representation for the solution of (||). 

To this aim it is convenient to rewrite the "first order" propagator S(x, y; A) 
in terms of a "second order" one 

S(x, y; A) = {i 1 v D v + m)A CT (x, y; A) , (19) 

A CT (x, y; A) being defined by 

(D^ + m 2 - \g a^F^W{x, y; A) = -5\x - y) , (20) 



G 



with a*" = J ; [7 M ,7 I/ ]. 

After replacing ( |l9| ) in (||), using an appropriate derivative it is. possible 
to take the differential operator out of the angle brackets and write fl] 

G ex (x 1 ,x 2 ;y 1 ,y 2 ) = -(ijidi^ + mi) (473 ^ + m 2 )H gl (x 1 ,x 2 ;yi,y2) ; (21) 
a similar expression can be given for G gI [x — y) , having set 

H gl (x!, x 2 ; yi, y 2 ) = --Tr c (U(x 2) xi)A^(xi,y-i;A)U(yi,y 2 )A 2 (x 2 ,y 2 ;-A)) , 
H&(x-y)=iTr c (U(y,x)A°(x,y;A)). (22) 

For the second order propagator we have the Feynmann-Schwinger repre- 
sentation 

/•OO fX />S -\ 

2 1 -2\ 



; p 00 px pS 

A CT (x,y;A) = — / ds Vzexp[-i / dr-(m 2 + i 
2 Jo Jy Jo 2 

5 s Pexp[z ff / drifA^z)], 
Jo 



(23) 



where the world-line connecting y to x is specified in the four-dimensional 
language by z^ = z^(r), in terms of an additional parameter r, and Sq = 

Texp - \ Jo drat*" ss Ji {z) and 8S^ = dz»Sz u - dz v 8z» (the functional 
derivative being defined through an arbitrary deformation, z — > z + 8z, of the 
world-line). 

Substituting (|23| ) into (^2|), we obtain covariant path integral represen- 
tations for H gl (xi, x 2 ; yi, y 2 ) and H gl (x — y). These representations involve 
again the gauge field only trough the Wilson correlator W qq for the four point 
function and a similar quantity W q for the two point function. W q is ob- 
tained replacing in (||) the loop T qq by a second loop made bw the world-line 
connecting y to x closed by the appropriate Schwinger string EfL 

For the Wilson loop correlator we can make the ansatz (|| -[ill) in the center 
of mass system even in the relativistic case, after taking C = (since in a 
relativistic treatment the perimeter term can be completely reabsorbed in a 
mass renormalization) . Notice that in such a case in general ( |ll| ) is only an 
approximation, however, if rewritten in an appropriate way, it seems to be a 



significant one. In fact, beside to have the correct-semi-relativistic limit, (11) 
becomes exact in important geometrical situations P] like those of a flat Wilson 
loop or of uniform rotatory motion. 

From the above path integral representations and by an appropriate re- 
currence method, an inhomogeneous Bethe-Salpeter equation and a Dyson- 
Schwinger equation can be derived for two other second order functions 
H(xi,x 2 ,yi,y 2 ) and H(x — y), which are simply related to H si (xi,x 2 ,yi,y 2 ) 
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and H gl (x — y), reduce to them in the limit of vanishing x\ — x 2 , yi — 2/2 or 
x — y and are completely equivalent for what concerns the determination of 
bound states, effective masses, etc. 

In the momentum space, the corresponding homogeneous BS-equation can 
be written (in a 4 x 4 matrix representation) 



t ■ / a ' u t t, l n k + u 1 _, k + u-, 

$p(fc) = — % I -r I ab (k-U,-P-\ ,-P ) 

V ' J (2tt) 4 V ' 2 2 ' 2 2 ' 

H X (\P+ k)<T a <S> P (u)a b H 2 ( -\p + k), (24) 

where &p(k) denotes an appropriate wave function, cr = 1, a = 0,/ij/ and the 
center of mass frame has to be understood (i.e. P = (ms,0), tub being the 
bound state mass). Similarly, the DS-equation can be written also 

f (fc) =I(^H k - i; ^ ^y*® ° b ■ w 

T(k) being the irreducible self-energy, defined by H(k) = Ho(k)+iHo(k)T(k)H(k). 

Notice that in principle (Q) and (25) are exact equations. However the 
kernels I a b are generated in the form of an expansion in a s and the string 
tension a. At the lowest order in both such constants, we have explicitly 

/ 0; o(Q;p,j/) - 16n^a sP a p' D aP (Q) + 

+4a J d 3 (e~^<\C\e( P() )e(p' Q ) J d\{p 2 Q p' 2 - [Xp' p T + (1 - \)p p' T } 2 }? 
d^e-^<e( P0 ) ^- C ^ 2 p' 

ICIvpo-Pt 

U tP a{Q;p,p') = ~^a s p a (S^ p Q a - 8^Q p )D afi {Q) + 

[ j3/- -iQ C CpPa ~ CcrPp I is 

i„„(Q\p,p) = tj».«0» - KQMS'Q, - €Q,)i>MQ) (26) 

where in the second and in the third equation (b = has to be understood. 
5 Reduction of the BS equation and spectrum. 

To find the qq spectrum, in principle one should solve first ( p5| ) and use the 
resulting propagator in (E4T) . In practice this turns out to be a difficult task and 
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one has to resort to the three dimensional equation which can be obtained from 
( p4| ) by the so called instantaneous approximation. This consists in replacing 
Hj(p) in j2~i| ) with the free quark propagator 2 ~* a and the kernel I a b with 

M st (k,k') obtained from I ab setting k = k' a = T " 2 Wl t w[ m } W2 t w * 

ao V ' 1 au o u u mi+m2 2 rai+m2 2 

with uij = y/"i 2 + k 2 and = .y/m| + k' 2 . 

The reduced equation takes the form of the eigenvalue equation for a 
squared mass operator, M 2 = Mq + £/ , with Mq = \J m\ + k 2 + a/ m 2 + k 2 
and 



(klf/lk') = ^J^^IOJ^^. (27) 
(27r) ,;i V 2wiW2 y 2w 1 w 2 

The quadratic form of the above equation obviously derives from the second 
order character of the formalism we have used. 

Alternatively, in more usual terms, one can look for the eigenvalue of the 
mass operator or center of mass hamiltonian i?cM = M = M + V with 
V defined by M V + VM + V 2 = U. Neglecting, consistently, the second 
order term V can be obtained from U simply by the kinematic replacement 

(w 1 +w 2 )(w' 1 +w' 2 ) i 



In ref. £§] we have evaluated the spectrum using both the operator M 2 
and Hqm, including the hyperfine terms in the potentials but omitting the 
spin-orbit ones. The numerical procedure we have followed was very sim- 
ple. It consisted in solving first the eigenvalue equation for the zero order 
hamiltonian Hq = Mq — + err by the Rayleigh-Ritz method, using the 
three-dimensional harmonic oscillator basis and diagonalising a 30 x 30 matrix. 
Then we have evaluated the quantities (Vv|-Hcm|Vv) and (?/v|M 2 |7/v) for the 
eigenfunctions ip v obtained in the first step. Since the quantity (V 2 ) was not 
completely negligible (ranging e. g. between few tens and 150 MeV for cc), 
we had to use slight different values of the adjustable parameters in the linear 
and in the quadratic formulation. We have kept the light quark masses fixed 
in both cases on typical current values, m u — m s — 10 MeV, m s — 200 MeV. 
Then for the calculations based on Hqm we have assumed m c = 1.40 GeV, 
rrib = 4.81 GeV, a B = 0.363 and a = 0.175 GeV 2 , taken essentially from the 
semi-relativistic fits. On the contrary for the calculations based on M 2 we 
have used a running coupling constant of the form 

" (ll-jN,)hg <28) 

cut at a maximum adjustable value a s (0) and with Nt = 4, A = 200 MeV. 
Furthermore we have chosen a s (0) = 0.35 and a = 0.2 GeV 2 in order to 
reproduce the correct J/^> — r\ c separation and the Regge trajectory slope; 
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Figure 1: Heavy-heavy and light- heavy quarkonium spectra. 
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Figure 2: mm spectrum and Regge trajectories. 
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m c = 1.394 GeV and mt — 4.763 GeV in order to obtain exactly the masses of 
J/* and T. 

In both cases, on the whole, the agreement with the data is good, not only 
for bottonium and charmonium (as in ordinary potential models) , but also for 
the light-light and light-heavy systems which are here essentially parameter 
free. Examples are reported in figs. 1 and 2. Circlet and dotted lines refer 
to the linear formulation, crosses and broken lines to the quadratic one. The 
only serious disagreement concerns the light pseudo-scalar mesons related to 
the chiral symmetry breaking problem. In fact in this case, as well known, a 
strict connection should exist between wave function in (|24|) and irreducible 
self energy in ( p5f ) and the use of the free propagator can not be correct. 
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